Reply to a Comment of M. Continentino on "Universally diverging Griineisen 
parameter and the magnetocaloric effect close to quantum critical points" 
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We show that the comment |cond-mat/0408217| by Continentino on our recent paper [PRL 91, 
066404 (2003)] reaches incorrect conclusions as the comment wrongly extrapolates from results valid 
close to a classical phase transition into the quantum critical regime. 



In his comment |l| and a previous preprint [3| Conti- 
nentino claims that our result for the quantum-critical 
Griineisen ratio F = ^ ^ j^-i/l''^) (where a and Cp arc 
thermal expansion and specific heat, respectively) should 
be replaced by F ^ y^i/'A^ where if) is the shift exponent. 
We explain below why this claim is incorrect above the 
upper critical dimension, i.e. d+ z > 4. Below the upper 
critical dimension, one has ip = i/z and no discrepancy 
arises. 

We will focus on the quantum critical regime, about 
which the comment raises the issue; the experiments re- 
ported in Ref. also concern this regime. Above the 
upper critical dimension {d + z > 4), a scaling analysis 
(see below) is more subtle due to the presence of some 
dangerously irrelevant variable(s), but a direct calcula- 
tion of F in the quantum critical regime is straightfor- 
ward. Therefore we reported in the results of explicit 
calculations for spin-densitv-wave (SDW) quantum crit- 
ical points (QCPs) [1,0, Hi J where the quartic couphng 
is dangerously irrelevant and '0 ^ vz. For instance, for 
d = 3 and z = 2 where the exponents ?/; = 2/3 and 
vz = 1, we find F ~ while Continentino would pre- 
dict F ^ T^"^/^. These calculations already show that 
the claim of Continentino is incorrect. 




control parameter r 



FIG. 1: Schematic phase diagram near a quantum-critical 
point. In the shaded area, classical scaling is expected. The 
temperature dependence in the quantum critical regime is 
measured along a path indicated by the arrow. The dashed 
lines, defined for positive r, describe the upper and lower 
crossover temperatures, Tcri oc (r/u)^ and Tct2 oc r'^' , re- 
spectively. 



It is straightforward to pin down why the rather gen- 
eral argument given in the comment fails. The reason 
is that both the Ehrcnfcst relation [ll and the scaling 
ansatz used in Ref. Q [Eq. (5) of Ref. |3; the ansatz is a 
special form of Eq. (1.30) in Ref. [^J are only valid close 
to the classical phase transition, i.e. in the shaded area 
shown in Fig. ^ but not in the quantum critical region. 
The part of the calculation of Ref. Q under discussion 
is, however, performed along the arrow depicted in Fig.^ 
which never enters the shaded region. 

The problem with Contincntino's approach can be 
traced back to the qualitatively different role which the 
quartic coupling u > of magnetic fluctuations plays in 
the classical critical and quantum critical regimes: it is 
a relevant coupling in the former but (dangerously) ir- 
relevant in the latter. This makes it incorrect to use the 
simple scaling form Eq. (5) of Ref. to describe all the 
regimes. Instead, one has to include explicitly the quartic 
coupling u in the scaling ansatz for the free energy 
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where r oc p—pc measures the distance from the QCP. We 
assume d+z > 4, d > 1/v and d ^ z to avoid logarithmic 
terms. In the quantum critical regimes QCl and QC2, 
\x\,y <C 1 and a; + y 3> y , one obtains (omitting 
all multiplying constant factors) f{x,y) « 1 -|- a; -|- y. 
(x and y appear separately in the Fermi liquid regime 
below QC2.) In contrast, the scaling ansatz 0, Q of 
Continentino / oc |a; + can only be valid within 

the Ginzburg region of the classical transition (shaded 

area in Fig. i.e. only for |x 4- y| <C y for d < 

4 (see, e.g., Ref. 0). The use of this formula outside 
of its raiige of applicability in region QCl (as done in 
Refs. 0, I3 and implicitly when comparing our results 
to the Ehrenfest relation in Ref. 0) implies that F oc 
"7^ ^ i, for r = 0. This result is wrong on two 
accounts. First, it entirely misses the leading term which 
corresponds to Gaussian critical fluctuations [cf. the "1" 
term in the expression for f{x,y) above]. Second, it is 
incorrect even for the sub-leading term for if) ^ vz as it 
depends explicitly on the exponent a which characterizes 
only the classical but not the quantum-critical transition 



2 



(in some cases, e.g. for d = 2 and Heisenberg symmetry, 
a classical transition does not even exist). It therefore 
also contradicts numerous results in the literature, e.g. 
Refs. Hi]. 

In his comment Q Continentino stated that Gaussian 
fluctuations should not be taken into account in the scal- 
ing analysis. This is incorrect, since Gaussian fluctua- 
tions are part of the singular contributions at the QCP 
and in fact the leading singular terms in the QCl and 
QC2 regions. As explained above, the dominance of the 
Gaussian fluctuations in QCl and QC2 is fully compat- 
ible (see e.g. Refs. 0, Hi) with both the existence of a 
Nccl line and ijj ^ vz due to the dangerously irrelevant 
variable u,. 

To summarize, the comment reaches incorrect con- 
clusions by extrapolating from results valid very close 
to the classical phase transition into the quantum crit- 
ical regime where they are not valid. Furthermore the 



scaling ansatz used in Refs. 0, Q incorrectly treats the 
dangerously irrelevant quartic coupling above the upper 
critical dimension. 
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